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; ABSTRACT 

o ' 

' Context. Thanks to the Venus Express Mission, new data on the properties of Venus could be obtained in particular concerning 
its rotation. 

Aims. In view of these upcoming results, the purpose of this paper is to determine and compare the major physical processes 
influencing the rotation of Venus, and more particularly the angular rotation rate. 
M— I , Methods. Applying models already used for the Earth, the effect of the triaxiality of a rigid Venus on its period of rotation 
(-H 1 are computed. Then the variations of Venus rotation caused by the elasticity, the atmosphere and the core of the planet are 
rS i' evaluated. 

I , Results. Although the largest irregularities of the rotation rate of the Earth at short time scales are caused by its atmosphere 
^ ' and elastic deformations, we show that the Venus ones are dominated by the tidal torque exerted by the Sun on its solid 
I body. Indeed, as Venus has a slow rotation, these effects have a large amplitude of 2 minutes of time (mn). These variations of 
^ ■ the rotation rate are larger than the one induced by atmospheric wind variations that can reach 25 — 50 seconds of time (s), 
depending on the simulation used. The variations due to the core effects which vary with its size between 3 and 20s are smaller. 
Compared to these effects, the influence of the elastic deformation cause by the zonal tidal potential is negligible. 
' ' Conclusions. As the variations of the rotation of Venus reported here are of the order 3mn peak to peak, they should influence 
' past, present and future observations providing further constraints on the planet internal structure and atmosphere. 

O ' 

' Key words. Venus rotation 
^. ■ 

1. Introduction and why its spin axis has a small obliquity (Goldstein 
T— I ' 1964; Carpenter 1964; Goldreich and Peale 1970; Lago and 
^ ; The study of the hregularities of the rotation of a planet Cazenave 1979; Dobrovoskis 1980; Yoder 1995; Correia 
. provides astronomers and geophysicists with physical con- ^nd Laskar 2001, 2003), others studied the possible res- 
traints on the models describing this planet and allows a onance between the Earth and Venus (Gold and Soter, 
better understanding of its global properties. Ahhough igjg. giUg 2005; Bazso et al., 2010). But few studies have 
Venus is the planet sharing the most similarities with the ^gen made comparing the major physical processes influ- 
Earth in terms of size and density, some of its characteris- g^cing the rotation of Venus at short time scale. Karatekin 
tics are poorly understood like its atmospheric winds and al. (2010) have shown that the variations of the rotation 
its superrotation. With its thick atmosphere, the determi- yenus due to the atmosphere should be approximatively 
nation of the period of rotation of Venus has been and stih iq seconds for the 117 days time span characterizing the 
is a very challenging task. Today thanks to Venus Express planet's solar day, but they neglected the main effect of 
(the first spacecraft orbiting Venus since the Magellan impact of gravitational torques on the solid body of 
Mission in 1994) and Earth-based radar measurements, planet. Here, applying different models already used 
the period of rotation of Venus can be revised and its for the Earth, the irregularities of the rotation of Venus are 
variations evaluated. So here for the first time a complete computed. After clarifying the link between the length of 
theoritical study of the variations of the rotation of Venus ^ay (LOD) and the rotation rate of the planet, the effects 
on a short time scale is presented as well as its implications of ^he triaxiality of a rigid Venus on its period of rotation 
concerning the observations. are evaluated. This model is based on Kinoshita's theory 
Many studies on the rotation of Venus have already (1977) and uses the Andoyer variables (Andoyer, 1923). As 
been made. Several authors studied this rotation on a long Venus has a very slow rotation (-243.020 d), these effects 
time scale to understand why Venus spin is retrograde 
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have a large amplitude (2 mn peak to peak) and could be 
observable as this is shown in Section [S] 

As we know, the variation of the period of rotation of 
the Earth of the order of 0.6 milliseconds (ms) for the sea- 
sonal component, are principaly due to the tidal zonal 
potential, the atmosphere and the oceans (Munk and 
MacDonald, 1960; Lambeck, 1980; Barnes et al., 1983). 
Hereafter quantifying the effect of the zonal component of 
the solar tides on Venus (Section|3]), the effects of its atmo- 
sphere (Section [5]) and the impact of its internal structure 
(Section ini) are presented. 

To conclude, a comparison of the properties and am- 
plitudes of these various processes is analysed (Section[7]). 
Their observability and their implications concerning the 
properties of Venus are also discussed. We find that, on 
the opposite of what is often assumed in the literature, 
the atmosphere should not be the most important effect 
leading to variations of the rotation rate of Venus. Indeed, 
the tidal torque exerted by the Sun on the solid body of 
the planet should be larger by more than a factor of two. 



2. Rotation rate and LOD definition 

A day is defined as the time between two consecutive cross- 
ings of a reference meridian by a reference point or body. 
While these periods vary with respect to time, we can use 
their average values for the purpose. For the Earth, the 
solar day (24 h) and the sidereal day (23.56 h) are defined 
respectively when the Sun and a star are taken as ref- 
erence. They are very close because the Earth's period of 
rotation is far smaller than its period of revolution around 
the Sun. By contrast, as Venus has a slow rotation, its so- 
lar and sideral days are quite different. Indeed, the canon- 
ical value taken for its rotational period (sideral day) is 
243. 02d (Konopliv et al, 1999), whereas its solar day varies 
around 117d. In the following, to be consistent with the 
physical inputs in the global circulation model (GCM) 
simulations, we will fix the value of the solar venusian day 
at 117d. 

There are mainly two methods to measure the rotation pe- 
riod, and its variations. The radar Doppler measurements 
give direct access to the instantaneous rotation rate uj{t) 
that we can translate in an instantaneous period of rota- 
tion LOD{t) by 

LOD{t) = ^ (1) 
w(i) 

Because we will be interested in variation around a mean 
value, we further define Aa;(t) = bj(t) —ZJ and ALOD{t) — 
LOD{t) - LOD where LOD = 243.02d and aJ = At 
first order ^{t) = -^^{t). Thanks to infrared im- 
ages of the surface of the planet, we can also measure the 
longitude of a reference point with respect to a given 
reference system at several epochs where (pi = io{t). As 
we will discuss in Section [71 as cj(i) has periodic varia- 
tions, fitting 4>i by a linear function could lead to a wrong 
estimation of w. 



3. Effect of the solid potential on the rotation 
rate 

A first and simplified model to describe the variation of 
the rotation of Venus is to consider that the atmosphere, 
mantle and core of Venus are rotating as a unique solid 
body. This model enables us to use Kinoshita's theory 
(1977) with the Andoyer Variables (1923) as those de- 
scribing the rotation of the Earth. The rotation of Venus 
is then described by three action variables (G, L, H) and 
their conjugate variables {g, l,h). G represents the ampli- 
tude of the angular momentum and L, H respectively its 
projections on the figure axis and on the inertial axis (axis 
of the reference plane) such as 



L — G cos J and H = G cos / 



(2) 



where / and J correspond respectively to the angles be- 
tween the angular momentum axis and the inertial axis 
and between the angular momentum axis and the figure 
axis. Here the figure axis and the inertial axis coincide re- 
spectively with the axis of the largest moment of inertia 
and the axis of the orbit of Venus at J2000.0. The angle 
J is yet unknown. As shown in appendix A, for plausi- 
ble value of this angle, its impact on the rotation is weak 
with respect to the other effects taken into account in this 
article. For sake of clarity, J is set to in the follow- 
ing. In this coordinate system, the Hamiltonian related to 
the rotational motion of Venus is (Cottereau and Souchay, 
2009): 

K = Fo+ E + E' + U. (3) 

Fq is the Hamiltonian for the free rotational motion de- 
fined by 



^" - 2^^ 



cos^ I 
B 



){G^-L^ 



2 G 



(4) 



where j4,i3,C are the principal moments of inertia of 
Venus. E^ E' are respectively the components related to 
the motion of the orbit of Venus which is caused by plane- 
tary perturbations (Kinoshita, 1977) and to the choice of 
the "departure point" as reference point (Cottereau and 
Souchay, 2009). \J is the disturbing potential of the Sun, 
considered as a point mass, and is given at first order by : 



GM' 



.2G - A-B. 



P2(sin(5) 



,A-B 



\PI(svwS) cos2q! 



(5) 



where Q is the gravitational constant, M' is the mass of 
de Sun, r is the distance between its barycenter and the 
barycenter of Venus, a and b are the planetocentric lon- 
gitude and latitude of the Sun with respect to the mean 
equator of Venus and a meridian origin (not to be con- 
fused with the usual equatorial coordinates defined with 
respect to the Earth). The P™ ^^^^ the classical Legendre 
functions given by: 



(-1)"'(1 - x2)- d"+"(a;2 - 1)" 



2"n! 



(6) 
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The haniiltonian equations are : 



dt 



{L,G,H) 



dK 



d{l,g,h) 



,h) 



dK 



d{L,G,H)- 



(7) 



(8) 



As the components wi and lu2 of the rotation of Venus 
are supposed at first approximation to be neghgible with 
respect to the component W3 along the figure axis, this 
yields : 



and 
d 



(9) 



(10) 



Splitting uj into its mean value uj and a variation Auj, uj 
is given by : 

l(CA.,.-f. (U, 

Notice here that Fq and E + E' do not contain g, so 
that the variations of the rotation of solid Venus are only 
caused by the tidal torque of the Sun included in the dis- 
turbing potential U. To make explicit the dependence on 
the variable g, we express ?7 as a function of the longitude 
A and the latitude /3 of the Sun with respect to the orbit 
of Venus at the date t with the transformations described 
by Kinoshita (1977) and based on the Jacobi polynomials 
such as : 



d{CAuj) 
di 



d QM' 
dg 



2C -A~B 



(-i(3cos2/-l) 



-sin^/cos2(A-ft.)) + 



e=±l 



4 

3 3 
-sin^/cos(2/ + 2g)+ ^ -(1 + ecos/)^ 



cos 2(A — h — el — eg) 



(12) 



where j3 is, by definition, equal to in this case. The first 
component at the right hand-side of Ea. lfT^ . does not de- 
pend on the variable g^ so that this yields: 



3 

- sin^ /cos(2/ 



2g) 



e=±l 



ecos/)^ 



cos 2(A — h — el — eg) 



\dt. 



(13) 



where has been replaced by • a and n being 
respectively the semi-major axis and the mean motion of 
Venus defined by the third Kepler law 



QM'. 



(14) 



DiflFerentiating with respect to g in the Eo. lfT^ . Aw is 
given by : 

A-B 



Auj = - / 3 



sin^/sin(2Z + 2g) 



+ U{l + ecosI)' 
sin 2(A — h — el ~ eg) 



dt. 



(15) 



Notice that we can write l+g ~ (j) (Kinoshita, 1977) where 
(j) is the angle between the axis which coincides with the 
smallest moment of inertia and a reference point arbitrarly 
chosen as the " departure point" on the orbit of Venus at 
the date t (Cottereau and Souchay, 2009). The position 
of this axis is given with respect to the origin meridian 
itself defined as the central peak in the crater Adriadne 
(Konopliv et al, 1999; Davies et al., 1992). In fact is 
the angle of proper rotation of the planet (</> « w). In 
the following, will be used instead of I + g. To solve 
analytically Ea.([T5|). the developments of (-)'^sin2(/) and 
(-)'^sin2(A — h ± (f>) are needed as a function of time 
through the variables M and Ls (respectively the mean 
anomaly and the mean longitude of Venus) taking the ec- 
centricity as a small parameter (Kinoshita, 1977). In a first 
approach, the orbit of Venus can be considered as circular 
(i.e : e=0, instead of : e=0.0068). This yields : 



Auj 



B-A 



sin^ / cos(20) 



1 

2</> 



--(1 + cos/)^ cos2(L5 



1 



1 



(1 - cos 7)^ cos2{Ls 



2Ls - 2(1) 
1 

2Ls + 20 



(16) 



Notice that this equation is similar to the equation given 
by Woolard (1953) in his theory of the Earth rotation, 
using a different formalism based on classical Euler angles. 





VenuH 


Period of revolution 


224. 70d (Simon et al, 1994) 


Obliquity 


2°. 6358 (Cottereau and Souchay, 2009) 


Mean period of rotation 


-243.020d (Konopliv et al., 1999) 


Triaxiality : ^j^j^ 


-1.647941 10"" (Konopliv et al, 1999; Yoder, 1995) 




0.33600 (Yoder, 1995) 



Table 1. Numerical values. and i?v are respectively 
the mass and the radius of Venus. 



Taking the numerical values of the Table [U the largest 
terms of the variation of — are 



Aijj 

UJ 



= 2.77 10"^cos(2Ls - 2(/)) 

-1-6.12 10"^ cos 20 

-1.99 lO""cos(2Ls-f-20) 



(17) 
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where 2Ls — 2(j), 2(p and 2Ls+2(p correspond to the leading 
terms with respective periods 58 d, 121.80 d and 1490d. 
As the two last terms scale as the square of the obliquity, 
they are significantly smaller than the term with argument 
2Ls — 2(f) (2.77 10^^). As Venus has a very slow rotation 
which appears in the scaling factor 3 ^ , the variations 
of the rotation rate due to the solid torque are larger than 
the Earth ones which correspond to amplitudes of ^ « 
10"^". FigH] shows the relative variations of the speed of 



3.x lO"' 
2.x 10^' 
1.x 



-1.x 10-' 
-2.x 10-' 
-3.x 10-' 




Venusian days 

Fig. 1. Variation of during a 4 venusian solar days 
time span (468d). One venusian solar day corresponds to 
11 7d. 



rotation of Venus due to the solid torque exerted by the 
Sun during a four venusian days time span (468d). The 
orientation of the bulge of Venus at i = is given with 
respect to the origin meridian. The mean orbital elements 
are given by Simon et al. (1994). The choice of the time 
span will allow us to compare in the following the different 
effects which act on the rotation of Venus. 

The physical meaning of the leading term with ar- 
gument 2Ls — 2(f) can be understood using a simple 
toy model. Consider a coplanar and circular orbit. At 
quadrupolar order, the gravitational potential created by 
the planet in its equatorial plane is equal to the one cre- 
ated by three point masses, one located at the center of the 
planet and of mass Aly — 2/x and two other ones symetri- 
cally located on the surface of the planet along the axis of 
smallest moment of inertia and of mass /i (see Figl2]) with 
M = being the planet mean radius. The torque 

exerted on the planet by the Sun can thus be computed 
using the forces exerted on these three points only. As 
shown in Fig|2l the sign of the resulting torque depends 
on the quadrant of the x, y plane the Sun is located in, 
and thus changes four times during a solar day, whose the 
corresponding argument is 2Ls — 2<p. 

To evaluate the influence of the eccentricity of Venus 
on its rotation rate, we show in FigI2] the residuals after 
substraction of 4^ when the eccentricity is taken into ac- 
count in the development of (f)^ sin 20 and {^)^ sin 2(A — 
h ± (f>) with respect to the simplified expression given by 
the Eq.([T7|. 







y 






Rv 






















equatorial plane 



Fig. 2. Diagram of the gravitational force exerted by the 
Sun in the frame corotating with Venus. The — and -I- 
signs give the sign of the net torque in each quadrant. 




Fig. 3. Influence of the eccentricity on the variation of the 
speed of rotation of Venus. 



We can observe that the eccentricity of Venus acts on 
^ with an amplitude of « 10~^ for the terms 2Ls + M - 
2(f) (46.34 d) and 2Ls - M-20 (78.86 d). These variations 
are smaller than the leading coefficient seen in Ea. (fT7|) by 
two orders and more important than the other coefficients 
with arguments 2(f) and 2(f) -I- 2Ls. 



From the relation — (t) 



ALOD I 



i^Qij -{t), the variations 
of the LOD(t) due to the torque of the Sun are 120 s (i.e 
2 mn or 0.0014 d). The implication of these results will be 
discussed in the following (Section[7]). Note that these vari- 
ations are larger than the uncertainties on the period of 
rotation measurements of Venus given by Magellan, that 
is to say 243.0200 ± 0.0002d in Konopliv et al. (1999), or 
243.0185 iO.OOOld in Davies et al. (1992) and close to the 
amount of their difference. 
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4. Effects of the elasticity on the rotation rate : 
deformation due to the zonal tidal potential 

For the Earth, an important variation of the speed of ro- 
tation is due to the zonal potential which causes temporal 
variations of the moment of inertia C . In this section, these 
zonal effects are evaluated for Venus, considered as a de- 
formable body. The zonal part of the potential exerted by 
the Sun on a point M at the surface of the planet is given 
by the classical formula 



V2,s{5m,5) 



(sin2 ,5. 



M 



(18) 



where S and d]\i represent respectively the planetocentric 
latitude of the Sun (the disturbing body) and of the point 
M and r the distance between the barycenter of the Sun 
and the Venus one. The corresponding bulge produced has 
a potential (Melchior, 1978): 



AF2,v(<5M,<5,r') = fc2V^2. 



„/3 ' 



(19) 



where ^2 and i?v are respectively the Love number and 
the radius of the planet. Differentiating the potential pro- 
duced by Venus at its surface (MacCullagh's formula), 
AV2,v can be expressed as a function of the principal mo- 
ments of inertia 



AF2,v(<5, r') = ^^{dC - dA){sm' S - i), 



(20) 



where we take dA = dB as the deformation is purely zonal. 
By identification we obtain: 



3 M' 
'4Mv 



dC - 1 
2 Mvi?^ ■ 



(21) 



Because this deformation does not induce a change in the 
volume of the isodensity surfaces in the planet to first or- 
der, dl = \dC+dA = iAC+AA = ( see Melchior, 1978, 
for a rigorous demonstration). The Euler's third equation 
is Clo = constant, to first order this yields — — 
Thus the variations of the angular rate of rotation of Venus 
are given by : 



Aw 

UJ 



AC 



V r 



3M' 



(sin" S 



r C 



(22) 



where Mm is the mass of Venus. Expressing sin" 5i as a 
function of / and A where /, A are given in Section [3] and 
represent the obliquity of Venus and the true longitude of 
the Sun, Ea. ip^ becomes 



Auj 



C My 



' sin / 



-(1 - cos2A) - 



(23) 



Using the developments of (^)3 and (^)3cos2A with re- 
spect to time from Cottereau and Souchay (2009), Fig|4] 
shows the relative variation of the speed of rotation of 
Venus due to the zonal potential. The numerical values 
of fc2 = 0.295 ± 0.066 and the ratio ^fm = 0.3360 are 
taken respectively from Konopliv and Yoder (1996) and 
Yoder (1995). The expressions found correspond to varia- 
tions of the LOD{t) of 0.019s peak to peak which is very 
small with respect to the contributions studied in Section 
131 So the zonal part of the potential on a non rigid Venus 
has a little infiuence on its very slow rotation, by contrast 
with the Earth for which these variations with annual and 
semi-annual component are not negligible (of the order of 
10^3 s) with respect to the 1 d rotation (Yoder et al., 1981; 
Souchay and Folgueira, 1998). Note that the semi-annual 
components are particularly small for Venus because its 
eccentricity and obliquity are much smaller than their re- 
spective value for the Earth. 



4.x 10-'° - 



2.x 10-'° - 



-2.x 10-'° 




Venusian days 

Fig. 4. due to the zonal potential during a 4 venusian 
days time span (468d). 



5. Atmospheric effects on the rotation of Venus 

As Venus has a denser atmosphere than the Earth, and 
as we know presently that the Earth atmosphere acts on 
its rotation in a significant manner (Lambeck, 1980), it 
will be interesting to study the corresponding effects on 
Venus. In this section, the core and the mantle of Venus are 
supposed to be rigidly coupled and its atmosphere rotates 
at a different rate. The variations of the speed of rotation 
of Venus due to its atmosphere are given by : 

where G, AG represent respectively the angular momen- 
tum of the rigid Venus and its variation due to the atmo- 
sphere AG = — AGatm- The angular momentum AGatm of 
the atmosphere can be split into two components : 

— The matter term Gm which is the product of lo with 
the inertia momentum of the atmosphere 

— The current term G„ which is due to the wind mo- 
tions with respect to the frame solidly rotating with 
the planet 
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From Lebonnois et al. (2010a) we have 



(1 + k'^)Gn + G„ 



9 

cos^ 9vgdhd9d(j). 



(25) 



where 0, (j), Ry stand respectively for the latitude, longi- 
tude and radius of Venus, we is the zonal wind and fcj 
is the load Love number of degree 2 (Karatekin et al., 
2010). Here to determine the variations of the angular 
momentum, two simulations made with the global circu- 
lation model (GCM) of the Laboratoire de Meteorologie 
Dynamique (LMD) (Lebonnois et al., 2010b) are used and 
compared. These simulations have been obtained with the 
LMD Venus General Circulation Model, using conditions 
similar to those presented in Lebonnois et al. (2010a), ex- 
cept for the boundary layer scheme. The first simulation 
{GCMl) was integrated from a zero wind state and is 
very close to the simulation published in Lebonnois et al. 
(2010a), though the winds in the deep atmosphere (0 to 
40 km altitude) are slightly higher, due to the updated 
boundary layer scheme. The second one {GGM2) was in- 
tegrated with initial winds in superrotation where the re- 
sulting winds in the deep atmosphere are close to observed 
values. The results presented in FiglS]were obtained after a 
200 days integration for GCMl and 130 days for GCM2. 
Here we calculate the seasonal variation of the angular 
velocity of Venus induced by the atmospheric winds to 
compare it with the variation discussed in section [3] The 
numerical values of the radius and the principal moment 
of inertia are taken in Cottereau and Souchay (2009). 



5.x 10-^ 


3 

-J -5.x 10-^ 
< 

-1.x 10-"* 
-1.5x10-'* 



Venusian days 

Fig. 5. due to the atmosphere during a 4 venusian 
days time span for the GCMl (red curve) and the GCM2 
(blue curve). 

Figure [5] shows the relative variation of the speed of 
rotation of Venus due to the atmosphere for a 4 venu- 
sian days time span (470 d) with the GCMl (red curve) 
and GCM2 (blue curve) . These variations have an ampli- 
tude of 1.26 10-s peak to peak for GCAfl and 2.44 IQ-^ 
respectively for GCM2. For the two models, the large am- 
plitude of 4^ of the order of 10^^ comes from the current 




term and depends on winds of Venus whereas the matter 
term acts on the rotation with an amplitude of the order 
of 10-1°. 

Using the definition of the LOD{t) given in Section [51 
the atmospheric contributions to Venus rotational speed 
correspond to peak to peak variations of the LOD{t) of 
27s with the GCMl and 51s with the GCM2. These values 
are consistent with the value of LOD{t) of the planet of 
7.9s given by Karatekin et al. (2009, 2010), who used the 
simulation presented in Lebonnois et al. (2010a). The dif- 
ferences are related to the amplitudes of the zonal winds 
in the region of maximum angular momentum (10-30 km 
of altitude), which vary between the different simulations. 
The most realistic values for these winds are obtained with 
the GCM2 simulation, where the winds are close to obser- 
vational data obtained from the Venera and Pioneer Venus 
missions (Schubert, 1983). We compared here the results 
given by GGMl and GCM2 to give an idea of the uncer- 
tainties in our present understanding of the atmospheric 
circulation and its modeling. In the following, only the 
values obtained with the model GCM2 will be compared 
to the other effects as it is thought to be the most realis- 
tic. Note that the atmospheric effects are nearly two times 
smaller than the solid effects described in Section [3l 

6. Core effects on the rotation of Venus 

The interior of Venus is probably liquid as inferred from 
the orbiting spacecraft data (Konopliv & Yoder 1996). 
The internal properties of Venus are expected to be like 
the Earth with a core radius around 3120 km (Yoder 
1995), but with a noteworthy difference because there is 
no dynamo effect on Venus (Nimmo 2002). In this Section, 
we investigate the impact of such a core on the rotational 
motion of Venus and especially on the variation of its 
LODit). For that purpose, we numerically integrate the 
rotational motion of Venus by taking into account the in- 
ertial pressure torque. The equations governing the rota- 
tional motion of two-layer Venus are the angular momen- 
tum balance for the whole body 



-ijj AH = V 



and for the core 



d< 



- A Jfc = 



(26) 



(27) 



(Moritz and Mueller 1987) where H, are respectively 
the angular momentum of Venus and of the core (see also 
e.g. Rambaux et al 2007). The vector F is the gravita- 
tional torque acting on Venus. Here, we assume that the 
core has a simple motion has suggested by Poincare in 
1910 and we neglect the core-mantle friction arising at 
the core-mantle boundary. Then, the rotational motion of 
each layers of Venus is integrated simultaneously with the 
gravitational torque due to the Sun acting on the triaxial 
figure of Venus. The orbital ephemerides DE421 (Folkner 
et al., 2008) are use for the purpose. 
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First we double-checked the good agreement between 
the numerical and the analytical solutions given by 
Ea. (fT7|) for a simplified model of one layer solid rigid 
Venus case. The very small differences obtained (at the 
level of a relative 10~^°) may result from the use of differ- 
ent ephemerides (DE421 for the numerical approach and 
VSOP87 for the analytical one). Then we applied the pro- 
cedure for the two layers case. As the moments of inertia 
Ic of the core are not yet constrained by measurements, we 
used internal models (Yoder, 1995) for which Ic is taken 
from 0.01 to 0.05, according to the size of the core. Notice 
that the flattening of the core is scaled to the flatten- 
ing of the mantle by the assumption that the distribu- 
tion of mass anomalies is the same. Table [5] shows the 
new amplitudes of 4^ for the 4 main oscillations with 
arguments 2Ls - 24^,2Ls ~ 2(t> + M, 2Ls - 24> - M, 20 
and with corresponding periods 58.37d, 46.34d, 78.86d, 
121. 51d. The presence of the fluid core allows a differential 
rotation of the mantle and of the interior of the planet. As 
a consequence in first approximation the mantle is decou- 
pled from the interior and presents amplitude of libration 
larger than in the solid case. The amplitude of the libra- 
tion increased with the size of the core. At maximum, 
for Ic — 0.05, we obtain an increase in the oscillation of 
2Ls — 2(j) of 17% corresponding to a variation of -|-20.4s. 



Period / 


58.37d 


46.34d 


78.86d 


121.51d 


Ic 


2Ls - 20 


2Ls - 20 + M 


2Ls - 20 - Af 


20 


0.01 


2.859979 


0.053456 


0.013032 


0.006328 


0.02 


2.950486 


0.055147 


0.013445 


0.006528 


0.03 


3.046908 


0.056950 


0.013884 


0.006741 


0.04 


3.149846 


0.058874 


0.014353 


0.006968 


0.05 


3.259981 


0.060932 


0.014855 


0.007211 




2.774866 


0.052244 


0.012726 


0.006120 



Table 2. Resulting amplitude for Venus with a fluid 
core (expressed in unit of 10"^). The last line is computed 
for a model without a fluid core. 



Figure ini shows the relative variation of the speed of ro- 
tation of Venus due to the core during a 4 venusian days 
time span with Ic — 0.05 (red curve) and Ic = 0.01 (blue 
curve) . These variations have peak to peak amplitudes re- 
spectively of 1.7 10~^ and 9.7 10^^. Using the definition of 
the LOD{t) given in Section [51 the core contributions to 
Venus rotational speed correspond to peak to peak vari- 
ations of the LOD{t) between 3.6s and 20.4s. Although 
the core effects on the rotation of Venus increase with 
the core size, they are smaller than the solid and atmo- 
spheric effects. Indeed, when Ic = 0.05 the core effects 
are respectively nearly two and six times smaller than the 
atmospheric and solid ones. 

7. Comparison and implication for observations 

The variations of the rotation of Venus presented in this 
paper are quasi-periodic and mainly due to three kinds of 
effects: solid, atmospheric and core. Compared to them, 
the zonal potential has a negligible influence. Finally Lo{t) 
can be written in the form: 

LL){t) = LJ + Auj 

= w+ [^y^(as,i + ac,i) cosjuJit + p^) 



+ 2J «aj- cos{ujjt + Pj) 



(28) 



where Wi, pi and as^i,ac,i (see Table are the frequen- 
cies, the phases and the corresponding amplitudes of the 
variation of rotation of Venus due to the solid and the 
core, and ojj, aaj (see Table S]) those due to the atmo- 
sphere. The atmospheric coefficients as well as their pe- 
riods have been obtained from a fast fourier transform 
(FFT) where pj are the phases. Note that the power spec- 
trum of the atmospheric variation is complex and can vary 
significantly from one model to another. As a consequence, 
only the most important frequencies are shown for the 
GCM2 model in Table IH 



t 

3 

-2.x 10-' 




Fig. 6. % 

span with 
curve) . 



Venusian days 

due to the core during a 4 venusian days time 
a /c = 0.05 (red curve) and a /c = 0.01 (blue 



Period Pi 


Argument 


solide as.i 




Core ac,i 


27V 




cos UJit 




cos CJit 


58d 


2L, - 20 


2.7710"" 


8.5110" 


< a, < 4.85 10"^ 


46.34d 


2Ls + M - 20 


5.2410"** 


1.01 10 


"" < a, < 8.510"" 


78.86d 


2Ls - M - 20 


1.2710"** 


3.110" 


1" <ae <2.110"'' 


121.80d 


20 


6.12 10"'' 


2.1 10" 


1° <a, < 1.110"« 



Table 3. Variation of ^ due to the solid and core effects. 



Comparing the amplitudes given in Tables [3] and |4j we 
see that the most important effect on the rotation rate 
of Venus is due to the solid potential exerted by the Sun 
on its rigid body. If all effects are taken into account, the 
variations of the LOD{t) can reach 3mn which could be 
observable in the future. 

As most past studies used infrared imaging of the sur- 
face to measure the evolution in time of the longitude of 
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Period Pj 


aa,j 


Pj 


2n 


COS U!jt + Pj 


117d 


4.17 10-'' 


3.11 


266d 


3.2210"^ 


1.74 


5.02d 


1.74 10"'^ 


2.39 


6.09d 


1.41 10"'^ 


3.14 


40.6d 


1.28 10"'^ 


0.0061 


79. 5d 


1.21 10"'^ 


1.64 



Table 4. Variations of ^ due to the atmosphere effects 
modeled by GCM2. 



reference points 0i at Venus surface, let us consider how 
this variable behave in our model. As 0i « uj{t) we have 



L(i) 



^1,0 



- Ljt - 



Aujdt. 



(29) 



FigEl shows the variations A0i during four venusian days 
time span, caused by the combined effects of the solid, the 
core and the atmosphere modeled by the GCM2. 




Venusian days 

Fig. 7. Variations of A0i x Ry in meter during a four 
venusian days time span. 

Like the variations of the rotation rate, these variations 
are periodic with a peak to peak amplitude reaching 12 m 
(time differences are converted in distance at the surface) . 
In precedent studies that measured the mean rotation rate 
of Venus, these variations where neglected. This implied 
fitting the measurements of the phase angle with a 
line of constant slope, this slope giving the mean rotation 
rate. Because the variations of the rotation rate discussed 
above are not negligible, we show in the following that 
this approach can yield large errors on the derived value 
of w, especially for a short interval of time. Similarly to 
Laskar and Simon (1988) let us consider the error made 
on the mean rotation rate when fitting the signal given by 
Eg. ip^ . keeping only the most important frequency for 
simplicity, by a function of the type : 



over a time span [Iq ,to + T]. For a least square fitting 
procedure, the residual is given by : 

pto+T 

-D = / {(j)i^Q + ult + As'mujit - {ujobst + (f>i.obs))'^dt 

Jto 

(31) 

where A and wi = ||rd/d correspond respectively to 
the larger amplitude of Eq. (|29p and its corresponding fre- 
quency. Minimizing this residual yields 



6A 



Tuji 



+2 



sin(a;iio) — sincji(T + to) 
2A 



cos(cjito) + cosa;i(T + ^o) 



(32) 



= ^ho + 7fr^ T{2T + 3to)iUiCOs{iUito) 
+T{T + 3to)LUi coswi (T + to) 
+3(T + 2io) [sin(cjito) - sin wi(T + to) 



(33) 



We can see that uiobs depends on both the time of the 
first observation (to) (i.e phase) and interval (T) between 
observations. If the phase of the effect is unknown, the 
max error made on the mean rotation rate is given by : 

ALOZJobs = max |wo6s(io, r) - (34) 
to 

FiglH] shows this maximum error ALODobs in seconds as 
a function of the interval T between observations. 




200 300 
T(Days) 



(30) 



Fig. 8. Maximum error of ALODobs{T) in seconds given 
by Ea.(IM|. The dotted line is the uncertainty given by 
Magellan (Davies, 1992) on the rotation of Venus. 



As we can see, even if the mean angular velocity is re- 
trieved for long baseline observations, the error yielded by 
the linear fit can be large for short duration observations. 
In addition, the use of such a simple model for pre- 
vents any measurement of the amplitude of the variations 
detailed in this paper. To measure the amplitude of these 
variations, modeling to and 0i with Eg. (1281) and Eq. 
respectively during the data reduction is necessary. 
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As the most important effect is due to the torque of 
the Sun on the Venus rigid body with a large amplitude 
on a 58d interval, it could be interesting to substract the 
measured signal by a fitted sinusoid of this frequency. The 
direct measurement of the amplitude of the sinusoid would 
give information on the triaxiality of Venus at 3 to 17% 
of error because of the core contribution. 

To disentangle atmospheric effects, a multi frequency 
analysis of the data will be necessary. Indeed after hav- 
ing determined the larger amplitude on 58d as explained 
previously, it could be substracted to the analysis. Then 
as the atmospheric winds (described by Eq.([55])) are the 
second most important effects on the rotation of Venus, 
the residuals obtained could constrain their strength. In 
parallel it should be interesting, also, to fit the signal ob- 
tained by a sinusoid on the period of 117d because only 
the atmosphere acts on the rotation with this periodicity. 
Note that this period is also present when using an al- 
ternative atmospheric model GCMl (Lebonnois, 2010a). 
So the corresponding amplitude could directly give indi- 
cations on the atmospheric winds. Fitting the signal with 
additional sinusoids with other periods given in Table S] 
could also increase the constraints. 

At last, as the core has not the same contribution 
on each period presented in Table [31 fitted the signal by 
Eg. ipS)) (or by Eq. ip^ ) could confirm the presence of a 
fluid core. Of course, such a detection is possible only if 
the precision of the measurements is of the same order of 
magnitude than the core effects and if the atmospheric ef- 
fects which add noise in the signal are better modelled in 
the future. 

Many values of the mean rotation of Venus have been 
estimated since 1975. FiglH] shows all these values as well 
as their error bars. The latter value of 243.023 ± O.OOld 
has been recently estimated from Venus express VIRTIS 
images (Mueller et al., submitted paper 2010). The dotted 
lines represent variations of 0.00197d around mean value 
of 243.020 ± 0.0002d (Konopliv et al., 1999) caused by 
solid, atmospheric (GCM2), and core (/c = 0.05) effects 
presented in this paper. Note that the value of Davies et 
al. (1992) set to 243.0185±0.0001dhas been recommended 
by the lAU (Seidelmann et al., 2002). 

As we can see from FigO 

— the variations of the period of rotation (0.00197d) pre- 
sented in this paper are consistent with most of the 
mean rotation periods measured so far. 

— if the true mean value of the period of rotation of Venus 
is close to the lAU value, the variations cannot explain 
the most recent value obtained by VIRTIS (Mueller 
et al., submitted paper 2010). Indeed, the difference 
between them of 7mn implies larger variations. 

The different values of the rotation of Venus since 1975 
could be explained by the variations of 0.00197d due to 
the solid, atmospheric and core effects. Despite the fact 
that the value of Davies et al. (1992) has been recom- 
mended by the lAU because of its small published error 
bars (±0.0001), its large difference with the recent VIRTIS 
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Fig. 9. Values of the period of rotation of Venus and their 
error bars given since 1975. The dotted lines represent 
variations of 0.00197d around the mean values of 243.020d 
caused by the effects presented in this paper. 

value (7mn) and the other measurements, is not in agree- 
ment with the variations presented in this paper. It seems 
difficult to explain these large variations, which would im- 
ply an increase of more than 50% of the effects discussed 
here, with the current models. So it would be interesting 
to compare in detail the different rotation rate measure- 
ment methods and to determine with a better accuracy 
the mean value of the rotation of Venus. 

8. Conclusion 

The purpose of this paper was to detemine and to compare 
the major physical processes influencing the angular speed 
of rotation of Venus. Applying different theories already 
used for the Earth, the variation of the rotation rate as 
well as of the LOD was evaluated. 

Applying the theory derived from Kinoshita (1977), 
the effect of the solid potential exerted by the Sun on a 
rigid Venus was computed. Considering in the first step 
that the orbit of Venus is circular, we found that the 
variations of the rotation rate have a large amplitude of 
2.7710"^ with argument 2Ls - 20 (58d). Taking into ac- 
count the eccentricity of the orbit adds periodic variations 
with a 10~* amplitude. On the opposite of the Earth, as 
Venus has a very slow rotation, the solid potential has a 
leading influence on the rotation rate which corresponds 
to peak to peak variations of the LOD of 120s 

Considering Venus as an elastic body, we then evalu- 
ated the impact of the zonal tidal potential of the Sun. 
These variations correspond to peak to peak variations on 
the LOD of 0.014s which are very small with respect to 
the contributions of the solid effect. 

Then we computed the effects of both the core and the 
atmosphere (modeled by two different simulations of the 
LMD global circulation model). According to our compu- 
tation the atmospheric and core contributions to Venus 
rotational speed correspond to peak to peak variations of 
the LOD of 25— 50s and 3.5— 20.4s respectively at different 
periods. Despite its thickness, the impact of the venusian 
atmosphere modeled by our most realistic simulation on 
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the rotation is 2.4 times smaller than the contribution of 
the solid torque exerted by the Sun. The variations of the 
LOD due to the core, which increase with its size, are still 
smaller. 

At last we have shown that the variations of lu and 
iiv^f*! which reach 3mn and 12m respectively, need to be 
taken into account in the reduction of the observations. 
Ignoring these variations could lead to an incorrect esti- 
mation of (D. With the steadily increasing precision of the 
measurements, carrying a frequency analysis of the data 
modeled by either Eqs. (P5|) or (1^ will hopefully enable 
to put physical constraints on the physical properties of 
Venus (triaxiality, atmosphere, core). 

To conclude, the variations shown in this paper, would 
explain different values of the mean rotation of Venus 
given since 1975. The difference of 7 mn between the lAU 
value (243.0185, Davies et al., 1992) and the VIRTIS one 
(243.023 ± 0.001, Mueller et al, submitted paper, 2010) 
implies larger variations, not found here, probably due to 
systematic errors in the measurements of the rotation. In 
view of these new results and the recent study of Mueller 
et al. (submitted paper, 2010), it would be interesting to 
compare the different methods of measurement of the ro- 
tation of Venus, to identify the error sources with a better 
accuracy and to revise the value of the mean rotation. 
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by Kinoshita (1977) without supposing J = such as: 



Pa (sin (5) = -{3cos^ J -1) 



-(3cos2/- l)P2(sin/3) 



sin2/sin(A - h)P^ {sin l3) 
sin 2/^2^ (sin cos 2(A - h) 



■ sin 2 J 



- sin 2/P2 (sin /?) cos g 



- i(l + ecos/)(-l + 2ecos/) 

e=±l 

Piisin 13) sin(A - h ~ eg) 

— -esin/(l + ecos/) 



e=±l 
32/ 



P2^(sin /3) cos(2A - 2h - eg) 
+ sin^ J 

- sin^ /P2(sin/3) cos2(? -f - ^ esin/ 



e=±l 



(1 + ecos /)P2^ (sin /3) sin(A -h- 2eg) - — 



^ (1 e cos /)^P|(sin /3) cos 2(A - /i - eg) 

(A.3) 



e=±l 



Appendix A: Influence of the angle J on the 
rotation rate 

In Section 121 we assumed that the angle J between the an- 
gular momentum axis and the figure axis of Venus can be 
neglected. Here we reject this hypothesis of coincidence of 
the poles and we evaluate the impact of J on the rotation. 
According to EqjS] the variations of Aw is given by : 



and 



— (CAuj) = — — . 
dV ' dg 



(A.l) 



where 



U = 



GM' 



,2C- A-B. 



P2(sin<5) 



]P2 (sin(5) cos2q! 



(A.2) 



As it was done in Section [3l we express P2(sin5) and 
P| (sin 6) cos 2a as functions of the longitude A and the 
latitude f3 of the Sun with the transformation described 



-(3cos2/- l)P2(sin/ 



P|(sin(5) cos 2q; = 3 sin^ J 

cos 2/ + 2 X! sin2/P2 (sin/3) sin(A ~h~ 2el) 
1 



e=±l 



■ sin^ I Pi {sin /S) cos 2{X - h - el) 
o 

- ^ psin J(l + pcos J) 
p=±i 

3 

- - sin2/P2(sin/3) cos{2 pi + g) 

- i(l + ecos/)(-l + 2ecos/) 
£=±1 ^ 

P2 (sin P) sin(A -h- 2pel - eg) 

— ^ iesin/(l + ecos/) 



e=±l 
32 



P2^(sin 13) cos(2A - 2h - 2pel - eg) 
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i(l + pcosJ)2 



p=±l 



3sin2 /P2(sin/3) 



cos{2l + 2pg) ^ esiii/(l + ecos/) 



e=±l 



P2 (sin /?) sin(A - /i - 2pel - 2eg) 
+ E ^(l + ecos/)2 



£=±1 
32 



P2 (sin ;3) cos 2(A -h- pel - eg) 



(A.4) 



Assuming that (3 ~ Q and removing the components which 
do not depend on the variable g this yields: 



d(CAw) 



d GM' /a\3 2C -A- B 



dt 




e=±l 

3 cos(2A -2h- eg) 



3.2.1 1 
- sm I — cos 2(7 

4 2 16 



^ (1 + e cos 7)^3 cos 2(A - ft, - eg) 

e=±l 

H — I psmJ[l + pcosJ) 

\p=±i 

[^sin2/icos(2p/ + 5) 
— -e sin /(I + e cos /) 

e=±l 

3cos(2A -2h- 2pel - eg) 

+ E J(l + PC0sJ)2[3sin2 4 
p=±i 

cos(2/ + 2pg)+ ^ -(l + ecos/)2 



£=±1 



3 cos 2(A — h — pel ~ eg) 



(A.5) 



where 2^ has been replaced by ■n?{-)^. 

The variations of the speed of rotation of Venus due to 
the solid torque exerted by the Sun are obtained by devel- 
oping Eq. (IA.5I) as functions of time through the variables 
M and Ls taking the eccentricity as a small parameter 
(Kinoshita, 1977). Fig lA.ll shows the residuals after sub- 
straction of ^ obtained by the Ea. (|A.5p when the angle 



J = is taken into account and e = with respect to the 
variations given by Ea. (fT7| in Section [3l We arbitrarily 
take the angle J = 0.5°. 




Venusian days 



Fig. A.l. Influence of the angle J = 0.5° on the variation 
of the speed of rotation of Venus. 



The influence of the angle J — 0.5° on the rotation 
with amplitudes of ~ 10^^ is smaller than the leading 
coefficients seen in Ea. (jl7p by three orders of magnitude 
and than the influence of the eccentricity of Venus by one 
order of magnitude. As J is probably much smaller than 
0.5° as it is the case for the Earth (1"), its influence on the 
rotation can be neglected with respect to the other effects 
taken into account in this paper. 
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